Librerías

library(readxl)
library(tidyverse)

Load data

Data <- read_excel("Data.xlsx")
x <- "Año"
y <- c("Lat", "Long", "Altitud", "TMAX", "TMIN")

1 All Species

tabla_final <- data.frame(
  "Variable"= character(),
  "Trend" = numeric(),
  "t" = numeric(),
  "p" = numeric(),
  "P95_max" = numeric(),
  "P95_min" = numeric(),
  "F" = numeric()
)
for (i in 1:5) {
  tabla <- data.frame(
    "Variable" = NA,
    "Trend" = NA,
    "t" = NA,
    "p" = NA,
    "P95_max" = NA,
    "P95_min" = NA,
    "F" = NA
  )
  tabla$Variable <- y[i]
  model_g <- lm(formula(paste(y[i], paste(x, collapse = "+"), sep = " ~ ")), data = Data)
  tabla$Trend <- model_g$coefficients[[2]]
  tabla$t <- summary(model_g)$coefficients[2, 3]
  tabla$p <- summary(model_g)$coefficients[2, 4]
  tabla$X95_max <-  confint(model_g, "Año", level = .95)[, 2]
  tabla$X95_min <-  confint(model_g, "Año", level = .95)[, 1]
  tabla$F <- summary(model_g)$fstatistic[1]
  tabla_final <- rbind(tabla_final, tabla)
}
DT::datatable(tabla_final,
              class = "nowrap",
              options = list(autoWidth = TRUE,
                             columnDefs = list(list(
                               width = "50%",
                               targets = 0))))

2 Indivudual species

spp <- unique(Data$Especie)

compare.coeff <- function(b_g,se_g,b_i,se_i){
  return((b_g-b_i)/sqrt(se_g^2+se_i^2))
}

tabla_ind <- data.frame(
  "Spp" = character(),
  "Variable"= character(),
  "Trend" = numeric(),
  "t" = numeric(),
  "p" = numeric(),
  "P95_max" = numeric(),
  "P95_min" = numeric(),
  "F" = numeric(),
  "Dif" = numeric()
)

for (n in 1:length(spp)) {
  ind <- filter(Data, Especie == spp[n])
  
  if (nrow(ind) > 10) {
    for (i in 1:5) {
      tryCatch({
        tabla <- data.frame(
          "Spp" = NA,
          "Variable" = NA,
          "Trend" = NA,
          "t" = NA,
          "p" = NA,
          "P95_max" = NA,
          "P95_min" = NA,
          "F" = NA,
          "Dif" = NA
        )
        #General
        model_g = lm(formula(paste(y[i], paste(
          x, collapse = "+"
        ), sep = " ~ ")), data = Data)
        
        tabla$Spp <- unique(ind[[1]])
        tabla$Variable <- y[i]
        model_i <-
          lm(formula(paste(y[i], paste(
            x, collapse = "+"
          ), sep = " ~ ")), data = ind)
        tabla$Trend <- model_i$coefficients[[2]]
        tabla$t <- summary(model_i)$coefficients[2, 3]
        tabla$p <- summary(model_i)$coefficients[2, 4]
        tabla$X95_max <-  confint(model_i, "Año", level = .95)[, 2]
        tabla$X95_min <-  confint(model_i, "Año", level = .95)[, 1]
        tabla$F <- summary(model_i)$fstatistic[1]
        
        
        
        b_g <- summary(model_g)$coefficients[2,1]
        se_g <- summary(model_g)$coefficients[2,2]
        b_i <- summary(model_i)$coefficients[2,1]
        se_i <- summary(model_i)$coefficients[2,2]
        
        
        tabla$Dif <- 2*pnorm(-abs(compare.coeff(b_g,se_g,b_i,se_i)))
        
        tabla_ind <- rbind(tabla_ind, tabla)
        
      }, error = function(e) {
        cat(
          paste0("WARNING: Specie ", ind[1, 1], " variable (", y[i], ") has"),
          conditionMessage(e),
          "\n"
        )
      })
    }
  } else{
    print(paste0("Data for ", ind[1, 1], " specie are insufficient"))
  }
}
## [1] "Data for Aegialia arenaria specie are insufficient"
## WARNING: Specie Agolius abdominalis pecoudi variable (Altitud) has 0 (non-NA) cases 
## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"
## WARNING: Specie Agolius bonvouloiri cantabricus variable (Altitud) has 0 (non-NA) cases 
## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"
## WARNING: Specie Ahermodontus marini variable (Altitud) has 0 (non-NA) cases 
## [1] "Data for Ammoecius brevis specie are insufficient"
## WARNING: Specie Ammoecius dentatus variable (Altitud) has subíndice fuera de  los límites 
## WARNING: Specie Brindalus porcicollis variable (Altitud) has subíndice fuera de  los límites 
## [1] "Data for Ceratophyus alloini specie are insufficient"
## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"
## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"
## [1] "Data for Chilothorax brancoi specie are insufficient"
## WARNING: Specie Chilothorax conspurcatus variable (Altitud) has 0 (non-NA) cases 
## [1] "Data for Chilothorax pictus specie are insufficient"
## [1] "Data for Erytus cognatus specie are insufficient"
## [1] "Data for Esymus sesquivittatus specie are insufficient"
## [1] "Data for Euoniticellus pallens specie are insufficient"
## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"
## [1] "Data for Jekelius bahilloi specie are insufficient"
## [1] "Data for Jekelius chersinus specie are insufficient"
## WARNING: Specie Jekelius hispanus variable (Altitud) has 0 (non-NA) cases 
## WARNING: Specie Jekelius intermedius variable (Altitud) has subíndice fuera de  los límites 
## [1] "Data for Limarus zenkeri specie are insufficient"
## [1] "Data for Liothorax plagiatus specie are insufficient"
## [1] "Data for Melinopterus pubescens specie are insufficient"
## WARNING: Specie Neagolius heydeni variable (Altitud) has subíndice fuera de  los límites 
## [1] "Data for Neagolius montanus specie are insufficient"
## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"
## [1] "Data for Onthophagus albarracinus specie are insufficient"
## [1] "Data for Onthophagus merdarius specie are insufficient"
## [1] "Data for Onthophagus nuchicomis specie are insufficient"
## [1] "Data for Onthophagus semicomis specie are insufficient"
## WARNING: Specie Oromus alpinus variable (Altitud) has 0 (non-NA) cases 
## [1] "Data for Oxyomus sylvestris specie are insufficient"
## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"
## [1] "Data for Plagiogonus arenarius specie are insufficient"
## WARNING: Specie Planolinus fasciatus variable (Altitud) has subíndice fuera de  los límites 
## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"
## WARNING: Specie Pseudacrossus sharpi variable (Altitud) has 0 (non-NA) cases 
## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"
## WARNING: Specie Rhyssemus parallelus variable (Altitud) has subíndice fuera de  los límites 
## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"
## WARNING: Specie Silphotrupes escorialensis opaculus variable (Altitud) has 0 (non-NA) cases 
## WARNING: Specie Silphotrupes orocantabricus variable (Altitud) has 0 (non-NA) cases 
## [1] "Data for Thorectes ferreri specie are insufficient"
## [1] "Data for Trypocopris vernalis specie are insufficient"

3 Variables

3.1 Latitude

3.1.1 Table

tabla_ind_lat <- filter(tabla_ind, Variable == "Lat")
DT::datatable(tabla_ind_lat,
              class = "nowrap",
              options = list(autoWidth = TRUE,
                             columnDefs = list(list(
                               width = "50%",
                               targets = 0))))

3.1.2 Graphs

for ( i in 1:length(spp)){
  ind <- filter(Data, Especie == spp[i])
  if (nrow(ind) > 10) {
  print(ggplot() + 
  geom_smooth(data= Data, aes(x = Año, y = Lat),col = "black", fill = "black", method = "lm") +
  geom_smooth(data= ind, aes(x = Año, y = Lat),col = "red", fill = "red", method = "lm")+
  ggtitle(paste0(spp[i]))+
  labs(x= "Year", y = "Latitude")+
  theme_minimal())
  }else{
  print(paste0("Data for ", ind[1, 1], " specie are insufficient")) 
  }
}

## [1] "Data for Aegialia arenaria specie are insufficient"

## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"

## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"

## [1] "Data for Ammoecius brevis specie are insufficient"

## [1] "Data for Ceratophyus alloini specie are insufficient"

## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"

## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"

## [1] "Data for Chilothorax brancoi specie are insufficient"

## [1] "Data for Chilothorax pictus specie are insufficient"

## [1] "Data for Erytus cognatus specie are insufficient"

## [1] "Data for Esymus sesquivittatus specie are insufficient"

## [1] "Data for Euoniticellus pallens specie are insufficient"

## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"

## [1] "Data for Jekelius bahilloi specie are insufficient"

## [1] "Data for Jekelius chersinus specie are insufficient"

## [1] "Data for Limarus zenkeri specie are insufficient"

## [1] "Data for Liothorax plagiatus specie are insufficient"

## [1] "Data for Melinopterus pubescens specie are insufficient"

## [1] "Data for Neagolius montanus specie are insufficient"

## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"

## [1] "Data for Onthophagus albarracinus specie are insufficient"

## [1] "Data for Onthophagus merdarius specie are insufficient"

## [1] "Data for Onthophagus nuchicomis specie are insufficient"

## [1] "Data for Onthophagus semicomis specie are insufficient"

## [1] "Data for Oxyomus sylvestris specie are insufficient"

## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"

## [1] "Data for Plagiogonus arenarius specie are insufficient"

## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"

## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"

## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"

## [1] "Data for Thorectes ferreri specie are insufficient"

## [1] "Data for Trypocopris vernalis specie are insufficient"

## Longitude ### Table

tabla_ind_lon <- filter(tabla_ind, Variable == "Long")
DT::datatable(tabla_ind_lon,
              class = "nowrap",
              options = list(autoWidth = TRUE,
                             columnDefs = list(list(
                               width = "50%",
                               targets = 0))))

3.1.3 Graphs

for ( i in 1:length(spp)){
  ind <- filter(Data, Especie == spp[i])
  if (nrow(ind) > 10) {
  print(ggplot() + 
  geom_smooth(data= Data, aes(x = Año, y = Long),col = "black", fill = "black", method = "lm") +
  geom_smooth(data= ind, aes(x = Año, y = Long),col = "red", fill = "red", method = "lm")+
  ggtitle(paste0(spp[i]))+
  labs(x= "Year", y = "Longitude")+
  theme_minimal())
  }else{
  print(paste0("Data for ", ind[1, 1], " specie are insufficient")) 
  }
}

## [1] "Data for Aegialia arenaria specie are insufficient"

## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"

## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"

## [1] "Data for Ammoecius brevis specie are insufficient"

## [1] "Data for Ceratophyus alloini specie are insufficient"

## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"

## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"

## [1] "Data for Chilothorax brancoi specie are insufficient"

## [1] "Data for Chilothorax pictus specie are insufficient"

## [1] "Data for Erytus cognatus specie are insufficient"

## [1] "Data for Esymus sesquivittatus specie are insufficient"

## [1] "Data for Euoniticellus pallens specie are insufficient"

## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"

## [1] "Data for Jekelius bahilloi specie are insufficient"

## [1] "Data for Jekelius chersinus specie are insufficient"

## [1] "Data for Limarus zenkeri specie are insufficient"

## [1] "Data for Liothorax plagiatus specie are insufficient"

## [1] "Data for Melinopterus pubescens specie are insufficient"

## [1] "Data for Neagolius montanus specie are insufficient"

## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"

## [1] "Data for Onthophagus albarracinus specie are insufficient"

## [1] "Data for Onthophagus merdarius specie are insufficient"

## [1] "Data for Onthophagus nuchicomis specie are insufficient"

## [1] "Data for Onthophagus semicomis specie are insufficient"

## [1] "Data for Oxyomus sylvestris specie are insufficient"

## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"

## [1] "Data for Plagiogonus arenarius specie are insufficient"

## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"

## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"

## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"

## [1] "Data for Thorectes ferreri specie are insufficient"

## [1] "Data for Trypocopris vernalis specie are insufficient"

## Elevation ### Table

tabla_ind_elev <- filter(tabla_ind, Variable == "Altitud")
DT::datatable(tabla_ind_elev,
              class = "nowrap",
              options = list(autoWidth = TRUE,
                             columnDefs = list(list(
                               width = "50%",
                               targets = 0))))

3.1.4 Graphs

for ( i in 1:length(spp)){
  ind <- filter(Data, Especie == spp[i])
  if (nrow(ind) > 10) {
  print(ggplot() + 
  geom_smooth(data= Data, aes(x = Año, y = Altitud),col = "black", fill = "black", method = "lm") +
  geom_smooth(data= ind, aes(x = Año, y = Altitud),col = "red", fill = "red", method = "lm")+
  ggtitle(paste0(spp[i]))+
  labs(x= "Year", y = "Elevation")+
  theme_minimal())
  }else{
  print(paste0("Data for ", ind[1, 1], " specie are insufficient")) 
  }
}

## [1] "Data for Aegialia arenaria specie are insufficient"

## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"

## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"

## [1] "Data for Ammoecius brevis specie are insufficient"

## [1] "Data for Ceratophyus alloini specie are insufficient"

## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"

## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"

## [1] "Data for Chilothorax brancoi specie are insufficient"

## [1] "Data for Chilothorax pictus specie are insufficient"

## [1] "Data for Erytus cognatus specie are insufficient"

## [1] "Data for Esymus sesquivittatus specie are insufficient"

## [1] "Data for Euoniticellus pallens specie are insufficient"

## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"

## [1] "Data for Jekelius bahilloi specie are insufficient"

## [1] "Data for Jekelius chersinus specie are insufficient"

## [1] "Data for Limarus zenkeri specie are insufficient"

## [1] "Data for Liothorax plagiatus specie are insufficient"

## [1] "Data for Melinopterus pubescens specie are insufficient"

## [1] "Data for Neagolius montanus specie are insufficient"

## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"

## [1] "Data for Onthophagus albarracinus specie are insufficient"

## [1] "Data for Onthophagus merdarius specie are insufficient"

## [1] "Data for Onthophagus nuchicomis specie are insufficient"

## [1] "Data for Onthophagus semicomis specie are insufficient"

## [1] "Data for Oxyomus sylvestris specie are insufficient"

## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"

## [1] "Data for Plagiogonus arenarius specie are insufficient"

## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"

## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"

## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"

## [1] "Data for Thorectes ferreri specie are insufficient"

## [1] "Data for Trypocopris vernalis specie are insufficient"

## Maximum temperature ### Table

tabla_ind_tmax <- filter(tabla_ind, Variable == "TMAX")
DT::datatable(tabla_ind_tmax,
              class = "nowrap",
              options = list(autoWidth = TRUE,
                             columnDefs = list(list(
                               width = "50%",
                               targets = 0))))

3.1.5 Graphs

for ( i in 1:length(spp)){
  ind <- filter(Data, Especie == spp[i])
  if (nrow(ind) > 10) {
  print(ggplot() + 
  geom_smooth(data= Data, aes(x = Año, y = TMAX),col = "black", fill = "black", method = "lm") +
  geom_smooth(data= ind, aes(x = Año, y = TMAX),col = "red", fill = "red", method = "lm")+
  ggtitle(paste0(spp[i]))+
  labs(x= "Year", y = "Maximum temperature")+
  theme_minimal())
  }else{
  print(paste0("Data for ", ind[1, 1], " specie are insufficient")) 
  }
}

## [1] "Data for Aegialia arenaria specie are insufficient"

## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"

## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"

## [1] "Data for Ammoecius brevis specie are insufficient"

## [1] "Data for Ceratophyus alloini specie are insufficient"

## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"

## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"

## [1] "Data for Chilothorax brancoi specie are insufficient"

## [1] "Data for Chilothorax pictus specie are insufficient"

## [1] "Data for Erytus cognatus specie are insufficient"

## [1] "Data for Esymus sesquivittatus specie are insufficient"

## [1] "Data for Euoniticellus pallens specie are insufficient"

## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"

## [1] "Data for Jekelius bahilloi specie are insufficient"

## [1] "Data for Jekelius chersinus specie are insufficient"

## [1] "Data for Limarus zenkeri specie are insufficient"

## [1] "Data for Liothorax plagiatus specie are insufficient"

## [1] "Data for Melinopterus pubescens specie are insufficient"

## [1] "Data for Neagolius montanus specie are insufficient"

## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"

## [1] "Data for Onthophagus albarracinus specie are insufficient"

## [1] "Data for Onthophagus merdarius specie are insufficient"

## [1] "Data for Onthophagus nuchicomis specie are insufficient"

## [1] "Data for Onthophagus semicomis specie are insufficient"

## [1] "Data for Oxyomus sylvestris specie are insufficient"

## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"

## [1] "Data for Plagiogonus arenarius specie are insufficient"

## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"

## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"

## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"

## [1] "Data for Thorectes ferreri specie are insufficient"

## [1] "Data for Trypocopris vernalis specie are insufficient"

## Minimum temperature ### Table

tabla_ind_tmin <- filter(tabla_ind, Variable == "TMIN")
DT::datatable(tabla_ind_lat,
              class = "nowrap",
              options = list(autoWidth = TRUE,
                             columnDefs = list(list(
                               width = "50%",
                               targets = 0))))

3.1.6 Graphs

for ( i in 1:length(spp)){
  ind <- filter(Data, Especie == spp[i])
  if (nrow(ind) > 10) {
  print(ggplot() + 
  geom_smooth(data= Data, aes(x = Año, y = TMIN),col = "black", fill = "black", method = "lm") +
  geom_smooth(data= ind, aes(x = Año, y = TMIN),col = "red", fill = "red", method = "lm")+
  ggtitle(paste0(spp[i]))+
  labs(x= "Year", y = "Minimum temperature")+
  theme_minimal())
  }else{
  print(paste0("Data for ", ind[1, 1], " specie are insufficient")) 
  }
}

## [1] "Data for Aegialia arenaria specie are insufficient"

## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"

## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"

## [1] "Data for Ammoecius brevis specie are insufficient"

## [1] "Data for Ceratophyus alloini specie are insufficient"

## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"

## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"

## [1] "Data for Chilothorax brancoi specie are insufficient"

## [1] "Data for Chilothorax pictus specie are insufficient"

## [1] "Data for Erytus cognatus specie are insufficient"

## [1] "Data for Esymus sesquivittatus specie are insufficient"

## [1] "Data for Euoniticellus pallens specie are insufficient"

## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"

## [1] "Data for Jekelius bahilloi specie are insufficient"

## [1] "Data for Jekelius chersinus specie are insufficient"

## [1] "Data for Limarus zenkeri specie are insufficient"

## [1] "Data for Liothorax plagiatus specie are insufficient"

## [1] "Data for Melinopterus pubescens specie are insufficient"

## [1] "Data for Neagolius montanus specie are insufficient"

## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"

## [1] "Data for Onthophagus albarracinus specie are insufficient"

## [1] "Data for Onthophagus merdarius specie are insufficient"

## [1] "Data for Onthophagus nuchicomis specie are insufficient"

## [1] "Data for Onthophagus semicomis specie are insufficient"

## [1] "Data for Oxyomus sylvestris specie are insufficient"

## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"

## [1] "Data for Plagiogonus arenarius specie are insufficient"

## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"

## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"

## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"

## [1] "Data for Thorectes ferreri specie are insufficient"

## [1] "Data for Trypocopris vernalis specie are insufficient"